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We present a Master Equation formulation based on a Markovian random walk model that ex- 
hibits sub-diffusion, classical diffusion and super-diffusion as a function of a single parameter. The 
non-classical diffusive behavior is generated by allowing for interactions between a population of 
walkers. At the macroscopic level, this gives rise to a nonlinear Fokker-Planck equation. The diffu- 
sive behavior is reflected not only in the mean-squared displacement {{r'^{t)) ^ with < 7 < 1.5) 
but also in the existence of self-similar scaling solutions of the Fokker-Planck equation. We give a 
physical interpretation of sub- and super-diffusion in terms of the attractive and repulsive interac- 
tions between the diffusing particles and we discuss analytically the limiting values of the exponent 
7. Simulations based on the Master Equation are shown to be in agreement with the analytical 
solutions of the nonlinear Fokker-Planck equation in all three diffusion regimes. 
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I. INTRODUCTION 

Diffusion is an ubiquitous phenomenon observed in physical, chemical, biological, social, algorithmic systems where 
"objects" (particles, molecules, cells, individuals, agents, ...) move in a seemingly random sequence of steps in such 
a way that their mean squared displacement increases linearly in time: (r^) ^ t where the proportionality factor is 
a constant which (apart from a numerical factor) is the diffusion coefficient. The resulting effect is a spread of the 
spatial distribution of the objects in the form of a Gaussian whose width grows with the square root of time. A 
microscopic mechanism yielding diffusive behavior at the macroscopic scale is based on the simple model of point 
particles undergoing random displacements on a one-dimensional lattice where they hop left or right from site to site 
at each tic of the clock: this is the random walk model for Brownian particles diffusion designed by Einstein who in 
1905 formulated the first microscopic derivation of the diffusion equation [1]. 

However there are many instances where the "objects" do not move freely: obstacles, time delays, interactions 
can modify their trajectories in such a way that the mean squared displacement deviates from the linear law and 
the Gaussian structure of the dispersion is deformed or replaced by a different distribution. So more generally, one 
observes (r^) ~ where for normal diffusion 7 = 1 while if 7 7^ 1 one talks about anomalous diffusion: when 7 < 1 
the process is said to be sub-diffusive and when 7 > 1 it is super-diffusive. As a result, there has been considerable 
interest in developing stochastic models capable of generating such behavior and the difficulty in constructing models 
depends on the requirements imposed. The most fundamental constraint is, of course, the need to reproduce a mean- 
squared displacement that exhibits power-law behavior as a function of time but this is a rather weak constraint that 
can be fulfilled in many ways. Perhaps the strongest demand that can be made is for the existence of self-similar 
solutions which is equivalent to demanding that all moments scale similarly, (r^") ^ t""'. This implies that the 
distribution has the form f{r,t) ~ tr"' l'^(l){r jV l"^) for some function (^{x) which is the case for classical diffusion. 
Another characteristic of classical diffusion is a type of universality wherein the details of the microscopic model can 
be changed without affecting the macroscopic process. In the Einstein model, the jumps can be of lengths greater 
than one with different probabilities including rests (jumps of length zero) without affecting the process (and, indeed, 
with the same diffusion constant). More generally, diverse microscopic dynamics can give rise to "diffusion" at the 
macroscopic level, but the underlying mechanisms may be quite different; for instance the distinction should be made 
between tracer motion where experimentally one follows trajectories of distinguishable particles seeded in an active 
medium [2, 3] and the motion of tagged particles which, while identical to the medium particles, are made observable 
by radioactive or fluorescent markers [4]. The latter case is referred to as molecular diffusion which is the dynamical 
process considered here. 
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FIG. 1. (Color online) The mean squared displacement as a function of time obtained from the solutions of the master equation 
for the cases of sub-, classical- and super-diffusion with a Gaussian initial condition. Other details of the simulation are 
described in Section IV. Only a single parameter varied is varied, a in Eq.(19), and the scaling exponents are predicted from 
the generalized diffusion equation, Eq.(6), to be 7 = -^^p^- The line segments are curves of the form (r^) = at'' and verify the 
diffusive scaling with the predicted exponents at long times. In this figure the mean-squared displacement is given in units of 
lattice-spacings, Sr, and the time in units of the fixed jump time, St. 

One class of diffusion models involves the use of memory or, equivalently, correlated noise since it is easy to see that 
if in two successive steps, a random walker preferentially follows the first step with one in the same direction (or the op- 
posite) direction the rate of diffusion can be dramatically altered. Non-Markovian dynamics is the mechanism behind 
the fractional Brownian motion [5, 6] and certain lattice models such as the elephant random walk [7, 8]. Similarly, 
the use of correlated noise in the generalized Langcvin equation [9] leads to the fractional Fokkcr-Planck equation [10] 
describing the phenomenology of anomalous diffusion in large ensembles [11] and for single trajectories [12]. 

One might wonder whether a Markovian random walk can give rise to anomalous diffusion. It is probably the case 
that any type of scaling of the mean-squared displacement is possible via the application of an external force (which 
is to say, jump probabilities or waiting times that depend on spatial position) but the existence of a field obviously 
corresponds to a very particular physical circumstance. Recently, we showed that sub-diffusive behavior could also be 
realized in a random walk model in which the walkers interact with one-another [13, 14] and which is indeed Markovian 
and contains no external force. Starting with a rather general ansatz for the way in which interactions affect the 
probability of jumps we derived from the master equation, via a multiscale expansion, a nonlinear Fokkcr-Planck 
equation, a generalization of the so-called porous medium equation. When this was combined with the requirement 
for diffusive-like self-similar solutions, it was shown that the certain limits of the jump probabilities (as explained in 
detail below) had to depend on the concentration of walkers in a specific power-law form. The resulting description 
of sub-diffusion was quite different from that given by non-Markovian models. For example, the fractional Fokker- 
Planck equation leads to a stretched-exponential distribution while the nonlinear model gives algebraic power-law 
distributions. This approach was also extended to the description of nonlinear reaction-diffusion systems [15]. 

In this paper, we further extend our previous work to include a description of super-diffusion. Indeed, we show 
that the same microscopic random walk model can exhibit sub-diffusion, classical diffusion and super-diffusion at the 
macroscopic scale as illustrated in Fig.l simply by varying a single parameter. In all cases, the behavior is "diffusive" 
in the strong sense of allowing for self-similar solutions. It will also be apparent that the model shows a kind of 
universality in which different microscopic dynamics yield the same macroscopic behavior. The ingredients needed 
to give this range of behavior are both nonlinearity and (spatial) non-locality. While our previous work showed that 
nonlinearity alone was sufficient to give rise to sub-diffusion, we find that spatial non-locality is critical to achieve 
super-diffusion within our framework. In some sense, this might be thought of as a complement or dual to the temporal 
non-locality of the non-Markovian approaches. 

In the next Section, we review our formalism and apply it to a particular form of nonlocal interaction between 
walkers. In Section III, we present direct numerical solutions of the master equation and compare to the predictions 
of the macroscopic nonlinear Fokkcr-Planck equation. Wc show that in the super-diffusive case, considerable care is 
needed in extracting the continuum limit from microscopic simulations but that the two are indeed in agreement. The 
paper concludes with a few comments on our results. 
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II. RANDOM WALK MODEL AND NONLINEAR FOKKER-PLANCK EQUATION 



A. The Master Equation 



Consider a walker moving on a one-dimensional lattice whose sites are labeled by integers j = —2 , —1, 0, 1, 2, ... . 
The formal expression of the microscopic dynamics describing the diffusive motion of Einstein's random walk model 
reads 

n*{r;t + l)=£_^-^'>n*{r + l;t)+£_^+^'>n*{r~l;t), (1) 

where the Boolean variable n*{r;t) = {0,1} denotes the occupation at time t of the site located at position r and 
^'^^^'> is a Boolean random variable controlling the particle jump between neighboring sites (^^~^-' + ^(+^) = 1); the 
superscript index (±1) indicates specifically that the jumps occur over one lattice distance. Extending the possible 
jump steps over the whole lattice, Eq.(l) becomes 

n*{r;t + l)= ^ ^U) n* {r - j;t) , (2) 

j = -oo 

with X]j=^oo ~ 1 (which condition prevents conflicting occupations at the arrival site). 

The mean field description follows by ensemble averaging Eq.(2) with (n*(r;<)) = f{r;t) and (^''■'•') = Pj, where j 
is the position index; using statistical independence of the ^'s and n*, we obtain 

-\-oo 

f{r;t + 5t)= Pj{r-jSr;t)f{r^jSr;t), (3) 

j=-oo 

where the distance, in lattice units, between neighboring sites is denoted by 5r and Pj{i) is the probability of a jump 
of j sites from site £ with J^j Pj{^) = 1- Setting j = ±1 and P±i = 1/2 in Eq.(3), we have the master equation for the 
usual random walk wherefrom Einstein derived the classical diffusion equation [1]. Equation (3) generalizes Einstein's 
master equation providing a formulation which in the hydrodynamic limit leads to the description of non-classical 
diffusion: when Pj{t) has a functional form pj F[f(r,t)] ~ Pj f°'~^ (where pj is a prescribed spatial distribution), 
Eq.(3) becomes the master equation leading to the description of nonlinear diffusion [13]. Equation (3) can also be 
rewritten as 

/(r, t^5t)^ f{r, t) = ^ (P,(r - j Sr) f{r - jSr, t) - P,{r) f{r, t)) , (4) 
j 

which expresses the rate of change of the particle distribution as the difference between the incoming and outgoing 
fluxes at location r. 



B. The Fokker- Planck equation 



We now introduce the further generalization that the jump probabilities Pj {£) take a functional form depending on 
both concentrations at the starting point and at the end point of the jump 



Pj^PjF[f{r^j6r,t)Jir,t)] 



with 



1, 



(5) 



where the probabilities pj are drawn from a prescribed distribution and the bounding condition < Pj < 1 imposes 
< F{x, y) < 1 as well as restrictions on the functional form of F{x, y). Under these conditions, multiscale expansion 
of the master equation was shown to give the generalized Fokker-Planck (or generalized Diffusion) equation [14] 
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with the notation 



and where 













dx 


f 


dx 






Sr 
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(7) 



2 St 




(8) 



are identified as the drift veioeity and the diffusion coefficient respectiveiy. Note that for F{x,y) ~ 1, Eq.(6) reduces 
to the classicai advection-difFusion equation. 

Since the function F {x, y) is defined in terms of the jump probabilities, it must be bounded, and so must satisfy 



0<xF{x,y) <1 and <yF{x,y) <l, V x, y G [0, 1] . 
It was shown in [14] that self-similar solutions are possible if and only if 

dxF{x,y) dxF(x,y) 



lim 



dx dy 

where the scaling exponent a is related to the diffusion exponent by 

2 

7 



a + l 



(9) 



(10) 



(11) 



Our previous work focused on local models for which F[x, y) = F{x) in which case the bounds given above (9) demand 
that a > 1 and therefore < 7 < 1, which is the correct formulation of sub-diffusion [16] but excludes the case of 
super-diffusion (7 > 1 i — > a < 1). 



III. SUB-DIFFUSION AND SUPER-DIFFUSION 
A. A model for sub- and super-diffusion 



Both types of anomalous diffusion (sub- and super-) can be described in a single formulation when the jump 
probabilities depend on the occupation probabilities at both the starting point and the end point of the jump by 
means of the ansatz 



F(x,?/; WsjWe) ^ UJsFs{x) +OJeFe{y) , 



(12) 



where and We are weighting factors relative to the functionals of the concentrations at the starting point and at the 
end point of the jump. Since an overall scale factor is unimportant, we can divide by and write a = uj^jixig without 
loss of generality. Finally, in order to assure that the constraint (9) is satisfied, we normalize to give the ansatz 



F{x,y;a) 



F{x)+aFiy) 
F{x)+F{y) ■ 



(13) 



The positivity arguments, x,y and the constraints (9) and F(x,y;a) > imply that < a < 1. We emphasize that 
this ansatz is motivated by simplicity and the necessary mathematical requirements. 

Considering the case that there is no drift (C = in Eq.(6)), and in order that the general formulation describe 
diffusion, we should have a scaling solution of the form / (r, t) = t^'^/'^cj) (r/P^^) ; it was shown in [14] that the scaling 
hypothesis demands that 



d d 
lim I -^xF {x, y; a) - -^xF {x, y; a, a) 
y^x \ c)x ay 



Kx° 



(14) 
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for some constant a. Using the model functional (13) in the l.h.s of (14) gives 

1 + a l — axF'{xa,Q) 



2 2 F{x]a,a) 

which is solved to yield 



Kx""' , (15) 



where B is an integration constant; reinserting (16) into (13), we find 

1+ 



B f 2K , , , 

F (x; a, a) = exp , (16) 

\1 — a a — I J 



2 K V -x" 



/ \ \y J \l — a Q — 1 , 

F(a;,y;a,a) = — • (17) 



1 + a exp 

+ (f)^exp(|i|i^^^i^) 
The natural limit: limct_i.i F{x, y;a) = 1 requires = i(l + a) so that 



T?f ^ l + aG{x,y;a,a) 

F[x,y;a,a) ^ -, 18) 

1 + G{x, y; a, a) 



with 



l + g 

G{x,y;a,a) ^ (-1 exp f ^ ^ ) ' 

It is clear that the limits a — I and a — >^ 1 give normal diffusion {F{x, y; a, a) = 1) and, for any finite value of a; > 0, 
we have < F{x, y;a,a) < 1 , V y G [0, 1] . Furthermore, one has that 

^rnF{x,y;a,a) = l^l'^ aSl '^^d lim F(x, y; a, a) = | J J (20) 

To provide some interpretation for this model, we note that 

d j^, \ n , ^ G{x,y;a,a) ( a;""^ - 1 ^ 

— F(x,y;a,«) = (1 + a)-— — 1, (21) 

ox (1 + G(a;, y; a, a)) \ x J 



d rpf \ , \ G{x,y;a,a) ( 1 - y" ^ v , . 

— F{x,y;a,a)^{l + a)——— ( ). (22) 

oy (1 + G(a;,y;a,a)) 



Since < y < 1, the signs of these derivatives are determined by the factors on the right: (x" ^ ~ l) and (l — y" ^) , 
and so depend on whether (a — 1) is positive or negative: 

d d 
a>l=^ —F{x,y]a,a) < < —F{x,y;a) 

d d 
a<l=^ Q^Fix, y; a, a) > > — F(a:, y; a) (23) 

In the first case, a > 1, the jump probability decreases with the concentration at the starting point and increases 
with the concentration at the arrival site; in other words the jump rate is reduced by putting more walkers at the 
origin and increased by putting more at the terminus of a jump: this is analogous to an attractive interaction. For 
a < 1, we have the reverse situation: the jump rate is increased by putting more walkers at the origin and decreased 
by putting more at the terminus, thus emulating a repulsive interaction. In the standard problem with all walkers 
at the origin at t = 0, the distribution decays monotonically away from the origin; thus, if the particles repel, the 
distribution expands faster (i.e. tends to a uniform distribution more quickly) whereas if they attract, then this 
attraction slows down the spread of the distribution. The physical interpretation is that attractive interactions give 
sub-diffusion and repulsive interactions give super-diffusion. 
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B. Explicit form of the scaling solutions 



The scaling solutions to the generalized diffusion equation, Eq.(6) with (18) and (19), are obtained following the 

1 {Srf 
1 St 



development given in ref.[14]. Defining D = ^^J2, and using the notation 



A 

we find: 

(i) for the super-diffusive case, a < 1, 



a- 1 



(l + a)(l + a) 



1 



where B is determined by normalization 

l=A-^/^B^ r (1+ 2/2) rfy 

J —oo 

and exists (is finite) provided 



a — 1 

whereas the mean-squared displacement (r^^ ~ I^ca given by 

(r^) = f r (1 + y2) '/^"-'^ dy 

J —OO 

and will only be finite if 



and the mean-squared displacement 



(r2) = f A-^'^B^ f y^ (1 - y^) '^^'^''^ dy 



(24) 



/(r,t)=i-'^/2i?V(-i)(i + ±^) , (25) 



"^^<0 ^ a>-l (26) 



3a — 1 „ 1 3 

-<0 =^ a>- =^ 7<-. (27) 

a ~ 1 6 I 

(ii) for the sub-diffusive case, a > 1, the distribution has finite support so that 

/ (., t) = t-7/25V(.-i) (l - I ^) e (l - I ^) , (28) 

with the normalization condition 



The continuum results therefore imply that an initial distribution of the form B' (1-I-Sa (r/w)^)^/'^" ^^8 (^1 + (''Z^^)^^ 

with B' determined by normalization and = =p for a ^ 1, will evolve self-similarly with mean-squared displacement 
increasing as V . 

IV. NUMERICAL SIMULATIONS 
A. Comparison between the master equation and the nonlinear Fokker-Planck equation 



Because the macroscopic description (the nonlinear Fokker-Planck equation) is derived from the microscopic dy- 
namics (i.e. the random walk model) by means of a well-defined but approximate multiscale expansion it is interesting 
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FIG. 2. (Color online) The spatial distribution as a function of time as determined by direct numerical solution of the master 
equation (symbols) and the analytic distribution, Eq.(28) for the case a = 0.5 corresponding to a scaling exponent of 7 = 4/3. 
All quantities are dimensionless (i.e. Sr = St = 1). In both cases, the initial distribution was the predicted scaling solution as 
described in the text: for the figure on the left, the initial width was ui = 10 while for the figure on the right it was w — 100. 



to compare the two to verify that they are indeed in agreement for the new ansatz proposed here. We therefore wish 
to compare the result of the direct, numeric, solution of the master equation to the continuum model, the generalized 
diffusion equation, Eq.(6), which predicts the existence of exact scaling solutions. (The explicit forms of which are 
given in the Appendix.) 

In the following, we use the elementary probabilities p±i — p for 1 < i < n for some values of p and n and have 
set a ~ 0. Following our previous work [14], we will use n = 2 (although in fact we have obtained similar results for 
one-step jumps, n = 1). We found that the simulations were increasingly sensitive to the value of p as a decreased: 
in fact, when this parameter is too large, and when a < 1 (super-diffusive case) the distribution determined from the 
master equation did not decay smoothly but, rather, showed an oscillatory structure in both space and time with a 
spatial period of several lattice sites. This is because in this case the walkers repel each other. So, if in an initial 
configuration the probability satisfies /„ > /„+i then so much more probability can get transferred from site n to 
site n + 1 than is transferred from n + 1 to n + 2 that one finds /„+i > /„. In the next time step, this causes a flow 
in the opposite direction. We therefore fixed on a value of p = 10""^ which is small enough to avoid any lattice-scale 
oscillatory behavior for a > 0.5 although we could use, e.g. p = 10^'' for a > 0.75. 

Even with this problem under control, we did not immediately obtain agreement between the lattice and continuum 
models. The left panel of Fig. 2 shows several snapshots of the distribution for the case a = 0.5 (7 = 4/3) starting 
with a distribution of width w = 10. While the two distributions are similar, there are significant differences between 
them. This was reflected in the mean-squared displacement from a run of 1 million updates of the master equation 
which appeared to behave as a power-law but with an exponent 7 of approximately 1.5 which is considerably larger 
than expected. Trying a wider initial distribution improved the agreement with the theoretical distribution (see Fig. 2; 
right panel) but still gave an exponent of about 1.45 for the mean-squared displacement. In order to test whether this 
is an equilibration effect, we used runs of 10 million time steps and extracted the mean squared displacement from 
each window of 1 million steps: the results are shown in Fig. 3. It is clear that several million time steps are required 
for the exponent of the mean-squared displacement to stabilize but even then, the final value is significantly above 
the theoretical value. Finally, in Fig. 4, we show the long-time value (from the final window of 1 million updates) as 
a function of the initial width of the distribution. Fitting to a function of the form 7 (w) = ^00 + ^ (fi^) gives for 
a = 0.5 a value of 700 = 1-34 and for a = 0.65 a value of 700 = 1-23 which compare well with the theoretical values 
of 1.33 and 1.21 respectively. 

Why should the result be so sensitive to the width of the distribution? Note that in order to obtain the continuum 
limit, we have to expand the distribution f (x + n) in terms of f{x). However, the form of the jump probability 
actually requires that we evaluate terms of the form 

/(-I) (x + n) - (x) = {a- l)r-2 (x) /' (x) n + ... 

The point is that with a < 1, and given that f{x) < 1, the coefficient /"~^ (x) can be very large so that the expansion 
is only reasonable if f'{x) is sufficiently small which, in turn, requires a wide enough distribution in units of lattice 
spacing. 



8 



1.6 



















- • 










• w=IO 
O w=32 
® w=100 




0=0.5 . 


8 


• 




® 




® 


2 O 
S ® (S 


8 i _ 


• 

o 




s 




8 


8 8 


8 


0=0.65 
8 S 
















0=0.75 


a=0.9 


1 . 1 . 1 . 1 



Time (units of 1 timesteps) 



10 



FIG. 3. (Color online) The scaling exponent, 7, as determined from fits of the numerical solution of the master equation. The 
fits were to the predicted functional form < >= ao(ai + 1)^ and were performed independently with successive windows of 
10® timesteps. The figure also shows results for different initial widths, w, (in units of lattice spacings) as described in the text. 
The lines are drawn at the predicted values 7 — jj—- 



y 1.3 




FIG. 4. (Color online) Gamma as a function of the inverse of the initial width (in lattice units). The lines are fits to 
7(10) = ao + ai (;^)"^. For a — 0.5, the fit gives 7(10 -> 00) = 1.339 and for a = 0.65, 7(10 00) — 1.221. The expected values 
of 7 based on the nonlinear Fokker-Planck equation are | = 1.333 and l| = 1.212, respectively. 



In summary, to get agreement with the continuum limit for the smallest values of a considered here - which is to 
say, for the strongest super-difFusivity - it is necessary to begin with very wide distributions and to perform sufficiently 
long runs. However, it is important to emphasize that these statements are quantified in terms of lattice spacings and 
time-steps: this has nothing to do with the physical width of the distribution and the physical duration of the process 
since both of these depend on the, as yet unspecified, Sr and St. As the distribution becomes wider and the times 
become longer, we can decrease both of these quantities so that the physical width of the distribution and time of 
the process remain constant (and as well maintaining the ratio ^ which determines the physical diffusion constant). 
Hence, the issues explored in this section pertain to reaching the continuum Umit and do not indicate any constraints 
on the physics of the model. 
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B. Sub- and Super-diffusion from the same model 



Finally, with an understanding of the issues involved in taking the continuum limit, we turn to a demonstration of 
sub- and super-diffusion from the same model. To do so, we have to address two additional issues. The first is that the 
choice a = is problematic for sub-diffusion. The reason is that if the distribution is zero for, say, r > tq then with 
a = the probability for a walker to jump from any site r < rg to a site r > tq is zero (see Eq.(20)). Even with an 
initial condition that, technically, has infinite support (such as a Gaussian), the distribution will be numerically zero 
outside some range so that in the diffusive case the maximum support is limited. To be able to perform arbitrarily 
long simulations over arbitrarily large lattices, we must have a > in the jump probabilities (Eq.(18). The second 
issue is purely a matter of the choice of problem: we wish to display both sub-diffusion and super-diffusion in the same 
system while only changing a single control parameter. Because the rates of the two processes are very different, the 
comparison cither necessitates very long runs for the sub-diffusive process or very restricted ranges of exponents for 
the super-diffusive process. Our compromise presented here was to take somewhat higher elementary probabilities, 
Pn = 10"'^, n G [—2, — f , 1, 2], a = 10^^ and 0.65 < a < 2 and beginning with an initial Gaussian with width fO lattice 
units. The results are shown in Fig.l where, after an initial transient period, the three regimes of diffusive behavior 
are clearly visible. 



V. CONCLUDING COMMENTS 



Beginning with a general formalism for describing the statistics of a population of interacting random walkers 
developed previously [f 3, f 4] we have constructed a particular class of interactions that can give rise to sub-diffusion, 
ordinary diffusion and super-diffusion with the adjustment of a single parameter. This flexibility depends on the 
introduction of non-local as well as non-linear interactions between the walkers. Not only docs the mean-squared 
displacement show the characteristic behavior of these three regimes, but the systems are diffusive in the strong sense 
of admitting self-similar solutions. Our model also demonstrates universality: not only does our particular ansatz give 
the same macroscopic behavior independent of the choice of the parameter a, but the sub-diffusive regime reproduces 
the same nonlinear Fokkcr-Planck equation as in our previous work that was based on a local interaction [f4]. 

We note that our calculations arc based on a particular ansatz for the concentration-dependence of the jump 
probabilities, Eci.(f3). The ansatz was motivated by simplicity and the requirements on the probabilities such as 
positivity and normalization. Otherwise, there is nothing special or unique about the choices made here: in particular, 
the fact that we allow for jumps of length 2 and the form of the ansatz itself. While we were able to characterize the 
effect of our model in terms of effective attraction or repulsion of the population of walkers, no attempt was made to 
relate it to microscopic interactions (e.g. particular force laws). Nevertheless, the underlying physics is simply that 
an effective repulsion between the walkers gives rise to super-diffusion while an attraction causes sub-diffusion. 

One interesting aspect of our results is that the range of scaling exponents that can be generated is restricted to 
< 7 < f .5. The upper limit is due to the nature of the self-similar solutions which are power-laws. As such, they 
only give finite result for a finite range of moments of the spatial variable. Normalization (existence of zeroth-order 
moment) demands that a > — 1 which is not restrictive in terms of the scaling exponent. However, existence of the 
mean-squared displacement itself demands that a > | thus giving the quoted limit on 7. In contrast, one typically 
expects that the upper limit for super-diffusion is 7 = 2 corresponding to ballistic motion. The difference is that 
ballistic motion is not-stochastic: an initial condition of a delta-function concentration of walkers at the origin would 
evolve, at best, as two delta-functions moving in opposite directions at constant speed with the only stochasticity 
occurring in the initial choice of direction. In contrast, the "free" limit of our model is the case of no interactions 
corresponding to ordinary diffusion. Moving away from ordinary diffusion requires turning on interactions and there 
is no reason to expect that this should ever lead to deterministic, non-stochastic behavior. As a matter of fact ballistic 
motion, although exhibiting mathematically a form corresponding to 7 = 2, originates from a physical mechanism 
(such as tracer dispersion [2]) different from molecular diffusion where free (ballistic) motion occurs at the microscopic 
scale only in the mean free path regime. 
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